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The results of the third phase of the Super-Kamiokande solar neutrino measurement are pre- 
sented and compared to the first and second phase results. With improved detector calibrations, 
a full detector simulation, and improved analysis methods, the systematic uncertainty on the total 
neutrino flux is estimated to be ±2.1%, which is about two thirds of the systematic uncertainty for 
the first phase of Super-Kamiokande. The observed 8 B solar flux in the 5.0 to 20 MeV total electron 
energy region is 2.32± 0.04 (stat.) ± 0.05 (sys.) xlO 6 cm _2 sec _1 under the assumption of pure 
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electron-flavor content, in agreement with previous measurements. A combined oscillation analysis 
is carried out using SK-I, II, and III data, and the results are also combined with the results of 
other solar neutrino experiments. The best-fit oscillation parameters are obtained to be sin 2 612 = 
0.30t°;g 2 (tan 2 9i2 = 0.42±%° t) and Arn^ = 6.2+19 x 10" 5 eV 2 . Combined with KamLAND results, 
the best-fit oscillation parameters are found to be sin 2 812 = 0.31 ± 0.01 (tan 2 612 = 0.44 ± 0.03) and 
Am 2 i = 7.6 ±0.2 x 10 _5 eV 2 . The 8 B neutrino flux obtained from global solar neutrino experiments 
is 5.3 ± 0.2(stat.+sys.)xl0 6 cm" 2 s _1 , while the 8 B flux becomes 5.1 ± 0.1(stat.+sys.) x 10 6 cm" 2 s _1 
by adding KamLAND result. In a three-flavor analysis combining all solar neutrino experiments, 
the upper limit of sin 2 #13 is 0.060 at 95% C.L.. After combination with KamLAND results, the 
upper limit of sin 2 #13 is found to be 0.059 at 95% C.L.. 

PACS numbers: 14.60.Pq 

I. INTRODUCTION II. SK-III PERFORMANCE 



The third phase of Super-Kamiokande (SK-III) began 
in October 2006 and ended in August 2008 when the elec- 
tronics were replaced. In SK-III, all 11129 PMTs have 
acrylic and Fiber Reinforced Plastic (FRP) PMT covers 
(blast shields) which were added at the start of SK-II in 
order to protect against propagating shock waves from 
PMT implosions. In the inner detector, the active pho- 
todetector coverage is 40% ( 40% in SK-I and 19% in 
SK-II). Thanks to detector improvements and superior 
analysis techniques, the SK-III's solar neutrino flux mea- 
surement is more precise than either SK-I's 2](SK before 
the accident) or SK-II's (SK with 46.5% of its PMTs) 
even with an exposure of only two years. In particular, 
the water purification system, event reconstruction and 
selection tools, as well as Monte Carlo detector simula- 
tion were improved. They will be explained in Chapter 
II and III in detail. 



In Chapter IV, the results of oscillation analyses are 
presented. By adding SK-III data, it was found that 
the energy spectrum and the time variation of solar neu- 
trinos obtained from our measurements favor only the 
large mixing angle solution (LMA) by constrai ning the 
8 B and hep neutrino flux to SNO NC flux [H Hg| and 
Standard Solar Model (SSM) prediction Q respectively. 
The first result of three-flavor neutrino oscillation analy- 
sis with the full SK data set will be shown as well as the 
two-flavor analysis result. In the last section of Chapter 
IV, the 8 B flux value obtained from the results of all so- 
lar neutrino experiments (global solar analysis) will be 
shown to compare with the prediction of the SSM. 
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A. Water system 

A major background for the solar neutrino observation 
at SK is the radioactivity from radon (Rn) from the U/Th 
decay chain in the water. The water in the detector is 
made from natural mine water using a very high perfor- 
mance water purification system. Even though the water 
is extremely pure there is still some Rn remaining. The 
Rn background events are very similar to solar events, 
so it is very difficult to remove them using only analysis 
tools. To reduce it, we have upgraded the system since 
the end of SK-I, including the addition of a new heat ex- 
changer and two reverse osmosis units during the SK-II 
and III periods. 

In addition, we investigated the water flow in the de- 
tector by intentionally injecting radon-enriched water. 
Tracing the resulting background events in time from 
this injected Rn, we found stagnation of water in the 
top and bottom of the detector volume, which increased 
the background. To counter this effect, we installed new 
pipes and changed the water flow. Previously, the wa- 
ter was supplied from the bottom of the inner detector 
(ID) and drained from the top of both the ID and outer 
detector (OD). Now, it is supplied from the ID bottom 
and drained from the top and bottom in OD and the top 
in ID with a total flow of 60 tons/hr, which is two times 
faster than before. This final setting has been in effect 
since August, 2007. As a result of these improvements, 
we have a central region with half of SK-I's background, 
enabling a lowering of the energy threshold. 

Note that the excessive background near the wall and 
bottom consisting of 7 rays due to the FRP cover also 
existed in SK-II. This background could not be reduced 
by improving the water system. 



B. Event Reconstruction 

1. Vertex 

The event vertex reconstruction for solar neutrino 
analysis performs a maximum likelihood fit to the tim- 
ing residuals of the Cherenkov signal as well as the dark 
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FIG. 1: Solid line shows the vertex resolution for SK-III as ' 1 1 1 1 1 1 1 — 

a function of the true total electron energy, while the dashed " "'"O 1( ^P 1500 

line shows that of SK-I. X POSltlOII (Cm) 



noise background for each testing vertex The vertex 
reconstruction method in SK-III was initially installed in 
SK-II and further improved over SK-II. It now has better 
resolution than SK-I. 

Figure [T] shows the vertex resolution for SK-III. The 
vertex resolution is defined as the distance from the true 
vertex position containing 68.3% of all reconstructed ver- 
tices. The vertex resolution in SK-I for 5 MeV electrons 
is 125 cm; here it is improved to 100 cm. 

A bias in vertex reconstruction is called the "vertex 
shift" . The vertex shift is defined as the vector from the 
averaged vertex of reconstructed events to that of the cor- 
responding simulated Monte Carlo (MC) events. Because 
the vertex shift results in events moving in or out of the 
fiducial volume, it represents one of the main systematic 
uncertainties for the solar neutrino flux measurement. 

The vertex shift is measured by placing a Ni-Cf gamma 
ray source Q at several positions inside the detector 
(Hereafter, the calibration using this Ni-Cf source is 
called "Ni calibration" or "Ni events"). The recon- 
structed data vertices at the fiducial volume edge were 
shifted more than 10 cm from the real source position in- 
ward toward the detector center, while those of the MC 
simulation were shifted less than 3 cm. It was found 
that this shift in data was due to an electronic effect of 
the relative hit timing within a wide range (~100 nsec). 
We measured the timing linearity by artificially shifting 
the common stop signal of individual TDCs for each hit 
channel. We found that a correction of -0.7% to the hit 
timing was required to restore linearity. After the correc- 
tion was applied, the vertex shift shortened significantly. 
Figure [5] shows the vertex shift in SK-III with timing cor- 
rection. The definition of x and z in Figure [5] (and other 
variables) is explained in Figure |31 



FIG. 2: Vertex shift of Ni calibration events. The origin of 
the arrows shows the true Ni source position and the direction 
indicates the averaged reconstructed position direction. The 
length of the arrow indicates the magnitude of the vertex 
shift. All vertex shifts are scaled by a factor of 20 to make 
them easier to see. 
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FIG. 3: Definition of detector coordinate system. 



2. Direction 

The direction reconstruction is based on the SK- 
I method: a likelihood function is used to compare 
Cherenkov ring patterns between data and MC simula- 
tion. An energy dependence is now included in the likeli- 
hood function for SK-III. The ring pattern distributions 
and their energy dependences are simulated for several 
energy ranges using electron MC simulation events. 

Figure @] shows the likelihood as a function of the re- 
constructed event energy and the opening angle between 
the reconstructed direction and the direction from the 
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FIG. 4: Likelihood value for reconstruction of event direction 
as a function of the reconstructed total electron energy and 
opening angle between the reconstructed direction and the 
direction from vertex to each hit PMT. 




10 15 
True total electron energy (MeV) 

FIG. 5: The solid line shows the angular resolution of SK- 
III as a function of the true total electron energy, while the 
dashed line shows that of SK-I. 



vertex to each hit PMT. Figure [5] shows the absolute an- 
gular resolution, which is defined as the angle of the cone 
around the true direction containing 68.3% of the recon- 
structed directions. For SK-III, the angular resolution is 
improved compared to SK-I by about 10% at 10 MeV 
and is close to the limit due to multiple Coulomb scat- 
tering of electrons. Note that the improvement of vertex 
reconstruction also contributes to the improvement of an- 
gular resolution, especially in the low energy region below 
6.5 MeV. 



3. Energy 

The reconstruction of event energy is similar to that 
for SK-I 0. The most important modifications with re- 
spect to SK-I are due to photo-cathode coverage and 
blast shields. Starting with the number of in-time hit 



PMTs (V50 coincident within 50 ns after the subtraction 
of time-of-flight (TOF) of Chcrcnkov photon from the 
reconstructed vertex to the hit PMT position), several 
corrections, described below, are made. The resulting 
effective hit sum N eff has less position dependence than 
V50. From N eff , we determine energy. This procedure is 
further outlined in Q, and is also explained in Q. 
The definition of N e ff is: 
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T" ^tail ^dark 
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where the explanations for the factors are as follows: 

Xi m . This factor estimates the effect of multiple photo- 
electrons in the i-th hit PMT. If an event occurs 
close to a detector wall and is directed towards the 
same wall, the Cherenkov cone does not have much 
distance to expand, and the observed number of 
hits is small. The correction Xi for this effect is 
defined as 



Xi = 



3.0, 



Xi < 1 
Xi = 1 



(2.2) 



where x { is the ratio of hit PMTs in a 3 x 3 PMT re- 
gion surrounding the i-th PMT to the total number 
of live PMTs in the same area. The — log(l — x.-,) 
term is the estimated number of photons per one 
PMT in that area and is determined from Poisson 
statistics. When Xi — 1, 3.0 is assigned to Xi. 

etaii- Some Cherenkov photons being scattered or re- 
flected arrive late at the PMT, and make late hits 
outside the 50 nsec time window. To correct the 
effect of the late hits, the term 



^tail 



VlOO ~ V 50 - Nalive X Rdark X 50 nSCC 



V 



(2.3) 



is added where iVioo is the maximum number of hits 
found by a 100 nsec sliding time window search. 

tdark'- This factor corrects for hits due to dark noise in 
the PMTs. 



^dark 



50nsec x N ahve x R d , 



ark 



v 50 



(2.4) 



where N„ 



is the number of all live inner detector 



No. 



No 



(ID) PMTs and Rdark is the measured dark rate for 
a given data taking period. 

-: This factor is for the time variation of the number 
'of dead PMTs. N aU is total number of PMTs and 
for SK-III it is 11129. 
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FIG. 6: Relation between N e ff and the true total electron 
energy (MeV) obtained from MC simulation. The dotted line 
in the upper figure shows the fitted polynomial function. The 
lower figure shows the deviation of the reconstructed energy 
from the polynomial function. 



FIG. 7: Energy resolution function obtained by electron MC 
simulation. Black points show one standard deviation for a 
Gaussian fit of the MC simulation divided by the true total 
electron energy, while the red (dashed) line shows a fit to 
a polynomial function. The black (dotted) line shows SK-I 
energy resolution. 



: This factor accounts for the direction-dependent 

photocathode coverage. S(0i,(f>i) is the effective 
photocathode area of the i-th hit PMT as viewed 
from the angles (9i,<pi) to take into account the 
shadowing of PMTs for glancing angles 0j. S 
is determined by MC simulation with the FRP 
PMT covers; the resulting N e ff is checked by an 
electron linear accelerator (LINAC) and an in-situ 
deuterium-tritium neutron generator (DT) calibra- 
tion data. 

exp(^|y): The water transparency is accounted for by 
this factor, where is the distance from the re- 
constructed vertex to the i-th hit PMT. X(t) is the 
measured water transparency for a given data tak- 
ing period. 

Gi(t): This factor adjusts the relative quantum effi- 
ciency of the PMTs. The differences in the quan- 
tum efficiency depend on the fabrication date of the 
PMTs. 

After determining N e ff, an event's energy in MeV can 
be calculated as a function of N e ff. The relation be- 
tween N e ff and MeV is obtained using mono-energetic 
electron MC simulated events as shown in Figure [5] The 
conversion function from N e f f to MeV is determined by 
fitting the relation with a fourth-order polynomial func- 
tion for the lower energy region (< 25MeV) and a first- 
order polynomial function for the higher energy region 
(>25MeV). 

The systematic uncertainty of the reconstructed energy 
is checked by TINAC and DT calibration. 

When calculating the energy for data events, the wa- 
ter transparency value, as determined by decay electrons 
from cosmic-ray muons, measured for six-day intervals, is 
used as an input parameter. For MC events, the change 

200 400 600 ° 



in water transparency and the relative quantum efficiency 
is simulated. 

The detector's energy resolution is well described by a 
Gaussian function. The energy resolution is described by 

a(E) = -0.123 + 0.376^ + 0.0349£, (2.5) 

in units of MeV (see Figure [7|). The SK-I resolution is 
cr = 0.2468 + 0.1492^ + 0.0690£ which is shown also 
in Figure [7J For low energy SK-III events, the energy 
resolution is improved by 5%, which is mainly due to the 
improved vertex reconstruction. 



C. Energy calibration 

As for SK-I, the primary instrument for energy cali- 
bration in SK-III is LINAC. A detailed discussion of the 
LINAC calibration methods can be found elsewhere [l2j . 
Single electrons are injected into the SK detector at vari- 
ous positions and at energies between 4.4 and 18.9 MeV. 
However, we could only take data with 4.4 MeV electrons 
at two positions, because the tuning of the electron beam 
is difficult for the lower energies. For this reason we did 
not include the 4.4 MeV data, and the lowest energy we 
included in this analysis is 4.7 MeV. The reconstructed 
energies of LINAC events are compared against those of 
the MC simulation to determine the energy scale. The 
absolute correction factor for PMT quantum efficiency 
was tuned to minimize any deviation between data and 
MC. The effect of the water transparency change on the 
energy scale was estimated as 0.22% by averaging over all 
energies and positions. The uncertainty of the electron 
beam energy, determined by a Ge detector measurement, 
is 0.21% (the same as for SK-I). 
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FIG. 8: The z dependence of the energy scale measured by 
LINAC calibration. The marker shows an average of (MC- 
DATA)/DATA over four energies at each data-taking point. 
The error bars show the RMS of (MC-DATA) /DATA for four 
energies. The filled circle markers are for R=4 m (x=-4 m, 
y=0 m) and the square markers are for r=12 m (x=-12 m, 
y=0 m). The open circle markers show the z dependence 
obtained by DT calibration. The dashed and dotted lines 
show ±1 and 0.5%, respectively. The edge of the fiducial 
volume is the same as the edge of the plot window (from - 
1610 cm to 1610 cm). 



In addition to the LINAC calibration, energy scale cal- 
ibration is done using 16 N produced with DT [13j . The 
generated 14.2 MeV neutrons exchange their charge with 
16 O. The produced 16 N decays into /3 and 7 with a half- 
life of 7.13 seconds. With a Q- value of 10.4 MeV, 16 N 
most probably decays into 16 O, an electron with maxi- 
mum energy 4.3 MeV and a 7 ray of energy 6.1 MeV. The 
peak value of the energy distribution is taken to evaluate 
the energy scale. DT data-taking is faster than LINAC 
data taking, so more positions can be checked. 

The position dependence of the energy scale system- 
atic uncertainty was estimated using only LINAC cali- 
bration data for SK-I and II, while for SK-III, in addi- 
tion DT calibration data are used to take into account 
the z-dependence and (^-dependence of the energy scale. 
The z-dependence is measured by LINAC, and the in- 
dependence is measured by DT. The LINAC can only 
take data at <j> = 180 degrees, whereas the DT generator 
can take data at five different positions at the same z 
position: center(r=4 m), = 0, 90, 180, and 270 degrees 
(r = 12m). 

Figure [5] shows the difference of the energy scale be- 
tween LINAC data and LINAC MC as a function of 
LINAC position. By averaging over all positions, the 
z-dependence is estimated to be 0.06% (the difference 
between r=4 m and r=12 m is also included here.) 

Figure [5] shows the difference of energy scale between 
DT data and DT MC. The vertical axis is normalized 
by the average of measurements at the same z positions. 
The mean value of the deviation from the average is taken 
as the (^-dependence of energy scale, which is 0.35%. A 



FIG. 9: The 4> dependence of the energy scale measured 
by DT calibration. The marker shows (MC-DATA) /DATA 
for each point and the error bars show statistical uncertainty. 
The circle, square, and triangle markers are for z=+12 m, m, 
and -12 m respectively. Zaverage is the average difference of 
energy scale between MC simulation and DT data obtained 
at the same z-position. The dashed and dotted lines show ±1 
and 0.5%, respectively. 
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FIG. 10: Energy scale difference between LINAC direction 
and other direction measured by DT calibration(LINAC di- 
rection range includes 76% of 7 MeV electrons). The solid 
and dotted lines show ±1 and 0.5%, respectively. 



resulting uncertainty of ± 0.35% for the overall position 
dependence is estimated. 

16 N decays allow directional studies of the energy scale 
which are not possible with the LINAC beam. The ob- 
served energy at several positions in the detector is com- 
pared with the MC-simulated energy and the difference 
is shown to agree with values obtained from LINAC data 
and MC. The 16 N energy scale difference between LINAC 
direction (downward direction) and the average of the 
other directions is estimated as ±0.25%, as shown in Fig- 
ure [TUl This difference is taken as the directional uncer- 
tainty of the energy scale. 

Finally, the energy scale uncertainty is calculated to be 
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Water transparency 
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LINAC energy sys. 


0.21 
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0.53 
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TABLE I: Systematic uncertainty of the energy scale. 

0.53% which is summarized in Table fl] This is slightly 
smaller than the SK-I estimated value of 0.64%. 

Energy resolutions of LINAC events are also compared 
for data and MC simulation. Figure [TT] shows the differ- 
ence of the energy resolution between data and MC as 
a function of the total electron energy. From Figure [TTJ 
±2.5% systematic uncertainty is assigned to energy cor- 
related systematic uncertainty for the spectrum measure- 
ment. 
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FIG. 11: Energy resolution difference between MC simula- 
tion and data as a function of energy, obtained by LINAC 
calibration. 



Quantitative estimates of trigger efficiencies are also 
obtained from 16 N data. The lowest hardware threshold 
setting has been in effect since April 2008. At this set- 
ting, the SK-III trigger achieved more than 99% efficiency 
at 4.5 MeV total electron energy. Before this time, the 
trigger efficiency was more than 99% at 5.0 MeV total 
electron energy. Figure [12] shows the trigger efficiencies 
of the lowest threshold period in SK-III. 



D. Optical calibration 

1. Light Propagation in Water 

For light propagation in water, a three-part model of 
light propagation consisting of absorption and two kinds 
of scattering is adopted for SK-III (as well as both SK-I 
and SK-II). In contrast to the previous phases, for SK-III 
the models are tuned using nitrogen/dye laser calibra- 
tion. We measured the attenuation length of scattering 



and absorption for four wavelengths (337, 365, 400 and 
420 nm) and tuned the water coefficients based on the 
measurement. The water coefficients are described by: 



a ray (\) = — [r 2 + — 



; (A) 



rn 
A 1 



UabsW = ^4 + "long (A, a) 
Based on the data of Feb. 2007 and 2008, we deter- 



mined 



and a. a long (A, a) was determined by a 



thirrFparty independent measurement introduced in SK- 
nl[|. Near the 400 nm region, we modified the relation 
using nitrogen/dye laser calibration. The crossing point 
of the SK and independent measurements [ijj is varied by 
another parameter (a). The results of the determination 
of these parameters are shown in Figure 1131 

The time dependence and position dependence of the 
water quality are described in the following. 

a. Time Dependence The time variations of the wa- 
ter coefficients and the water transparency are measured 
simultaneously by nitrogen/dye laser calibration and us- 
ing decay electrons from cosmic-ray muons. These mea- 
surements confirmed that the change of the water trans- 
parency is mainly caused by the change of the absorption 
coefficient. We obtained the relation of the absorption 
coefficient and the water transparency in SK-III using 
those data as well as water transparency as measured by 
decay electrons from cosmic-ray muons, as for SK-I. 

Figure [H] shows the time variation of the measured 
water transparency during SK-III and the stability of the 
peak energy of the decay electrons in SK-III as a function 
of time. The stability of the energy scale has 0.47% RMS. 

b. Position Dependence It was found that the PMT 
hit rates measured for a Ni-Cf gamma-ray source at the 
top region in the detector are systematically lower than 
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FIG. 13: Wavelength dependence of the water coefficients: 
scattering and absorption combined (black, solid), absorption 
(red, dashed), Rayleigh scattering (green, dotted) and Mie 
scattering (blue, dashed-dotted). The absorption coefficient is 
a function of water transparency and the water transparency 
in this graph is 139 m. 
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FIG. 15: Hit rate difference, (MC-DATA) /DATA, for barrel 
(upper), top (left lower), and bottom (right lower). Black 
(dashed) shows /3=0. Red (solid) shows the hit rate difference 
with tuned (3. This calibration data was taken in Feb. 2008, 
the same period as the LINAC calibration. 
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FIG. 14: The upper figure shows the time variation of the 
measured water transparency (weighted by the Cherenkov 
spectrum) during SK-III. The lower figure shows the stability 
of the peak energy of decay electrons in SK-III as a function 
of time. 



of the absorption parameter is considered, because the 
dominant contribution to the time variation of the water 
transparency is the absorption. The depth-dependence 
can be modeled as 



a abs {X){l+ /3 ■ z), 
a„fc,(A)(l-0 - 1200), 



for z > -1200cm 
for z < -1200cm 
(2.6) 

Below z=-1200 cm, the absorption coefficient is assumed 
to be uniform due to convection of the water. (3 is called 
the TBA parameter (unit is cm -1 ) and it varies in time, 
because TBA varies. The f3 for each period was tuned to 
minimize the TBA for that period. 

Figure [T5] shows a comparison of the hit-rate between 
data and MC simulation with the tuned (3 value, which is 
8.85xl0 -5 cm -1 , and without j3 correction, respectively. 
The hit-rate is defined as the number of hits in each PMT 
in units of the averaged number of hits for all PMTs 
during the calibration run. As shown by the figure, data 
and MC show better agreement with the (3 correction. 



2. Reflectivity of Black Sheet 



those for the bottom region by 3 ~ 5%. This rate dif- 
ference is denoted the "top-bottom asymmetry" (TBA). 
The MC simulation cannot reproduce the top-bottom 
asymmetry with a uniform attenuation length through- 
out the detector volume. In order to solve this prob- 
lem, a simple model of the light absorption is introduced 
to take into account the dependence of the attenuation 
length on depth. In this modeling, depth-dependence 



In the detector simulation, the reflectivity of the black 
sheet which covers the ID wall is calculated using the 
law of Fresnel, taking into account the polarization of 
Cherenkov photons. To better model the real situation, 
we measure the reflectivity for three incident angles us- 
ing calibration data taken from a movable light injector 
and a black sheet reflector. We found that a value of 
half of the SK-I reflectivity gives better agreement and 
a wavelength-dependent correction is newly applied for 
SK-III. 
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FIG. 16: Timing distribution of each hit PMT after tuning 
the MC simulation, for LINAC 5.1 MeV (x,z)=(-4 m,0 m) 
data. The blue histogram with open circles shows data and 
the red dashed histogram shows MC simulation. 



In addition, to better describe the diffuse and specu- 
lar reflection on the black sheet, we use two models for 
reflection: the Lambert model is used to describe diffuse 
reflection, and the Phong model is used for description 
of specular highlights 0. 

3. PMT and Electronics 

The PMT must be simulated as precisely as possible. 
Reflection and quantum efficiency are tuned using the ni- 
trogen/dye laser calibration data. In contrast to the de- 
scription used for SK-I, we put a wavelength dependence 
into the reflection and an incident angle dependence into 
the quantum efficiency. The position and width of the 
PMT after-pulses and the timing resolution were tuned 
using data from single electrons injected into the SK de- 
tector by LINAC. After these tunings, the timing distri- 
butions of LINAC data and MC simulation agree very 
well, as shown in Figure [TBI 

III. DATA ANALYSIS 

A. Data set 

As for SK-I and II, SK-III has two trigger levels for 
solar neutrino analysis: low-energy (LE) and super-low- 
energy (SLE) thresholds. SLE triggered events are fil- 
tered online to reduce the amount of data written to lim- 
ited storage space. Events reconstructed outside the fidu- 
cial volume are rejected. The data are reduced by a factor 
of approximately six. From Aug. 4th, 2006 to Jan. 24th, 
2007, only the LE trigger threshold was applied. The LE 
trigger is 100% efficient above 6.5 MeV and the livetime 
of this period is 121.7 days. The threshold was lowered 
and the first SLE data were taken with more than 99% 



efficiency at 5.0 MeV from Jan. 24th, 2007 to Apr. 17th, 
2008. This period is called the SLE1 period, and the 
livetime of this period is 331.5 days. From Apr. 17th, 
2008 to Aug. 18th, 2008, a lower SLE trigger threshold 
was applied and the SLE trigger had more than 98% ef- 
ficiency at 4.5 MeV. This period is called SLE2 period, 
and the livetime of this period is 94.8 days. 

Two neutrino samples are used for SK-III analysis. 
The first sample, for event energies between 6.5 and 20 
MeV, has a total livetime of 547.9 days. The second sam- 
ple, for event energies between 5.0 and 6.5 MeV, has a to- 
tal livetime of 298.2 days after rejecting high background 
periods caused by radioactive impurities accidentally in- 
jected into the detector. 

B. Event selection 

Basic explanations of each selection step are as follows: 

• Run selection 

In this selection, short runs (< 5 min.), runs with 
hardware and/or software problems, or calibration 
runs are rejected. 

• Cosmic-ray muon 

Cosmic-ray muon events are rejected by a total 
charge cut (< 2000 photo-electron (p.e.) in ID). 

• Electronic noise reduction 

Events due to electronic noise are rejected. 

• Fiducial volume cut 

Events which have reconstructed vertex position 
within 2 m of the ID wall are also rejected. 

• Spallation cut 

This cut is to reject events caused by cosmic-ray 
muons. 

• Event quality cut 

In this step, results of event reconstructions are 
tested. In particular, (a) quality of vertex and di- 
rection reconstruction, and (b) hit pattern, (c) re- 
sult of the second vertex reconstruction used for 
SK-I 2} are checked. Events produced by flasher 
PMTs are also rejected in this step. 

• External event cut 

This cut is to reject events induced by radioactivity 
from the PMTs or the detector wall structure. 

• Cosmogenic 16 N cut 

Events caused by decay of 16 N are rejected. The 
16 N is produced when cosmic ray [i is captured by 
16 O in water. 

• Small clustered hit cut and tighter fiducial volume 
cut 

These cuts reject events which have clustered hits 
(see|mC]and[niD|. 
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For more detailed descriptions of the reduction steps, 
please refer to @, || . 



C. Small clustered hit cut for the lowest energy 
region 

This cut is newly developed for SK-III to reduce the 
low energy background in regions near the edge in the 
detector. As described in the previous section, the target 
background is assumed to be triggered by a coincidence 
of dark hits and small clusters of hits due to radioactive 
sources in the FRP or the structure of the detector wall. 

We can separate this background from the solar neu- 
trino signal by searching for a small cluster in both space 
and time. A real neutrino signal at the edge region also 
has similar characteristics, but it causes a bigger cluster 
compared to the background events. Thus, the key is to 
evaluate the size of the hit cluster. 

In order to check the size of the hit cluster for an event, 
a new cut variable was created. This variable is a product 
of two factors. One factor is the radius of a circle which 
contains the clustered hits in a event, and the other is the 
number of clustered hits within a 20 nsec time window. 
Figure [17] shows the distribution of this cut variable for 
real events and for solar neutrino MC simulation events in 
a energy region 5.0 to 6.5 MeV. For this figure, a volume 
cut with r > 12 m and z > — 3 m is applied to study 
background events in the barrel, which are due to FRP 
cover radioactivity. 

As expected, the real data sample which contains 
mostly FRP events shows smaller values for this cut 
variable than simulated solar neutrino events. A cut 
value of 75 is selected which gives maximum significance 
(Signal / \J BG) for the solar neutrino signal. This cut is 
applied to events with energy < 6.5 MeV whose vertex 
position is in an edge region. The edge region is deter- 
mined for 5.0-5.5 MeV, and 5.5-6.5 MeV energy bins sep- 
arately, to optimize the significance of the solar neutrino 
signal. The criteria of the edge region are determined as 
follows: 



X} 

£ 0.08 



CO 
CD 



0.02 









I 1 I 


' ' 1 


, , , | , , , , 




w c 


A 


N 




J 4 
1 f 
I ? 




k 

l <& 

•. $ 

\ 




/ p 

J, J: , 


, , 1 





50 100 150 200 

Cut variable 



FIG. 17: Distribution of cut variable for the small clustered 
hit cut. Blue (filled circle) shows background sample events 
which are selected from r > 12 m and z > —3 m. Red (open 
circle) shows the solar neutrino MC simulation and the same 
volume cut is applied. For both data and MC events, event 
selection cuts up to the external event cut are applied. The 
reconstructed energies of these samples are in 5.0-6.5 MeV. 
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FIG. 18: Vertex distribution of 5.0-6.5 MeV energy region 
before and after the small cluster hits cut. The horizontal 
axis shows r 2 and the vertical axis shows the number of events. 
Events with z > —7 m are selected to show background events 
in the barrel. 



( r 2 > 180 m 2 for 5.5 <E< 6.5 MeV , , 

\ r 2 > 155 m 2 for 5.0 <E< 5.5 MeV { ' 

where r is defined in Figure[3] Figure IT51 shows the vertex 
distribution of the final data sample for the 5.0-6.5 MeV 
energy region before and after the small clustered hits 
cut. 



D. Tighter fiducial volume cut 

This cut is to reject the remaining background in the 
edge region. As shown in Figure 1191 background events 
in the bottom region still remain in the final data sam- 
ple. This non uniformity of the vertex distribution of the 



background distorts the angular distribution of the back- 
ground, which causes a large systematic uncertainty for 
the day-night asymmetry of solar neutrino flux. 

To set a tight fiducial volume, the significance as a 
function of detector radius is calculated. Based on the 
significance calculation, the final value of fiducial volume 
is obtained for each energy region. 

5.0 - 5.5 MeV : (r 2 < 180 m 2 and z > -7.5 m) = 13.3 kton 
5.5 — 20 MeV : no tight fiducial volume cut = 22.5 kton 

Note that in the future, we hope to remove this tighter 
fiducial volume cut by improving our signal extraction 
method, because solar neutrino events are preferentially 
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FIG. 19: Vertex distribution of final data sample in the 5.0- 
6.5 MeV energy region (22.5 kton fiducial volume). 
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FIG. 20: Number of events after each event reduction step as 
a function of energy. From the top, the markers show events 
after spallation cut, event quality cut, external event cut, 16 N 
cut, and small clustered hit cut plus tight fiducial volume cut. 
The dashed histogram shows the number of events in the SK-I 
final sample. 



rejected in the low energy region by this cut. 



E. Summary of event selection 

Figure (50] shows the energy spectrum after each step 
of the reduction and Figure [21] shows the remaining ef- 
ficiency of 8 B solar neutrino MC with respect to the re- 
constructed energy. The number of events after each re- 
duction step is summarized in Table |TT] (real data) and 
IIIII (solar neutrino MC simulation). While the event rate 
in the real data as a function of energy is the same as 
for SK-I, the cut efficiencies are improved by 10% in the 
final SK III data sample. 
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FIG. 21: Reduction efficiency for the solar neutrino MC sim- 
ulation. The definition of each histogram is the same as for 
Figure [201 



Reduction step 


5.0-6.5 MeV 6.5-20.0 MeV 


Spallation cut 


1861770 


114350 


Event quality cuts 


734843 


58187 


External event cut 


54820 


43146 


16 N cut 


54351 


39879 


Small cluster hits cut 


42916 




Tight fidv. cut 


24311 




Final 


24311 


39879 



TABLE II: Reduction results for data. 



F. Simulation of solar neutrinos 

The method to extract the energy spectrum and the 
flux of solar neutrinos in SK-III is based on that of SK-I 
and II. The 8 B and hep neutrino spectra are generated 
separately. The total 8 B and hep flux values are referred 
from the Standard Solar Model (SSM) 0. We use the 
spectrum and its uncertainty of 8 B and hep neutrinos 
calculated by Winter [l(| and Bahcall [ll| respectively. 
The event time of a solar neutrino event is simulated 
from the livetime of the SK-III full operation period, so 
that the expected zenith angle distribution of the solar 



Reduction step 


5.0-6.5 MeV 6.5-20.0 MeV 


Total 


100 


100 


Bad run cut 


86 


87 


Trigger condition 


72 


86 


Flasher events cut 


56 


67 


Spallation+ 16 N cut 


45 


53 


Event quality cuts 


45 


46 


External event cut 


28 


42 


Small cluster hits cut 


27 




Tight fidv cuts 


24 




Final 


24 


42 



TABLE III: Reduction results for solar neutrino MC events 
in %. For the spallation cut and 16 N cut, position-dependent 
dead time is considered. 
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neutrinos can be simulated correctly. 



G. Systematic uncertainty 

The following items are updated with respect to SK-I 
with the improved calibration, detector simulation and 
analysis tools described above. 

• Angular resolution 

In order to estimate the systematic uncertainty on 
the total flux due to the angular resolution differ- 
ence between data and MC simulation, a predicted 
energy spectrum is made by artificially shifting the 
reconstructed direction of the solar neutrino MC 
events. The shifted direction is calculated as the re- 
constructed direction ± the systematic uncertainty 
of angular resolution with respect to the generated 
direction of the recoil electron. After the solar angle 
fitting, the systematic uncertainty due to the uncer- 
tainty of angular resolution is estimated as ±0.67% 
on the total flux in 5.0-20 MeV region (±1.2% for 
SK-I). This method of estimation is the same as for 
SK-I, so the improvement is due to the new detec- 
tor simulation and angular reconstruction. 

• Vertex shift and resolution 

The vertex shift difference between data and MC 
is used to estimate the systematic uncertainty of 
the fiducial volume. We measure this difference us- 
ing Ni calibration rather than the LIN AC, because 
LINAC electrons always point downwards while so- 
lar neutrino recoil electrons go in all directions. 
Also, very few positions near the edge of the fidu- 
cial volume can be probed by the LINAC. As an 
example, Figure [22] shows the observed vertex shift 
difference as a function of energy at (15.2 m, -0.7 m, 
12 m), where the vertex shift in the x direction is 
larger than other positions in Figured There is no 
evidence that the vertex shift difference is energy 
dependent. 

To estimate the total flux uncertainty due to the 
vertex shift, the reconstructed vertex positions of 
the solar neutrino MC events are artificially shifted 
outward. Then, the fraction of events rejected by 
the fiducial volume cut due to the shift is estimated. 
Figure [53] shows the systematic uncertainty on the 
fiducial volume due to the vertex shift which is 
based on the results of the Ni calibration. Note 
that in SK-I, and II, the estimation was done by 
a calculation based on the size of the vertex shift, 
whereas in SK-III, MC simulations after event se- 
lection are used to take into account the cuts' ef- 
ficiencies at the edge of the fiducial volume. This 
new method gives a more accurate estimation. 

This study shows that 0.54% of the MC events are 
rejected after shifting the vertex position, and this 
fraction is set to the total flux systematic uncer- 
tainty due to the vertex shift (±1.3% for SK-I). 
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FIG. 22: The vertex shift as a function of energy. The mark- 
ers show the amount of vertex shift measured by Ni calibra- 
tion at x=15.2 m y=-0.7 m and z=12 m. For the vertical axis, 
negative values mean vertex shifts inward towards the center 
of the detector. 
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FIG. 23: The systematic uncertainty of the fiducial volume 
as a function of energy. 



The improvement with respect to SK I is mainly 
due to the reduction of the vertex shift. 

In Figure I231 the step in the fiducial volume un- 
certainty between 5.0 and 6.0 MeV is due to the 
tight fiducial volume cut for the lower energy region 
(< 5.5MeV). This relative difference of fiducial vol- 
ume systematic uncertainty is taken into account in 
the uncertainty of the energy spectrum shape. Be- 
tween 5.5 and 20 MeV, since vertices are assumed 
to be shifted in the same direction for all energies, 
the relative differences are 0.1%. Consequently, the 
systematic uncertainty of the spectrum shape is set 
to 0.1% in that energy region. In the 5.0-5.5 MeV 
region, the systematic error is set to 0.5% because 
the uncertainty of the fiducial volume is relatively 
0.5% larger than for the other region. 

The vertex resolution is also compared for data 
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FIG. 24: The difference of vertex resolution between MC sim- 
ulation and data as a function of energy obtained by LINAC 
calibration. The dashed line and the dotted line show 1 cm 
and 2 cm difference respectively. 



and MC simulation using LINAC events. Figure [24] 
shows the difference of vertex resolution as a func- 
tion of energy. This difference results in only a 
second order effect on the systematic uncertainty 
of the fiducial volume. 

• Event quality cuts 

The systematic uncertainties associated with the 
event quality cuts are estimated separately for (a) , 
(b) and (c) of UTTbI For (a) and (b), the sys- 
tematic uncertainties on the total flux are taken 
from the differences of the efficiencies between the 
LINAC data and MC simulation which are ±0.4% 
and ±0.25%, respectively. 

For (c), the second vertex cut, the estimation is 
done in the same way as for the first vertex fitter, 
resulting in a 0.45% contribution to the total flux 
uncertainty. 

The combined uncertainty is 0.65% which is about 
one third of SK-Fs. The improvement is due to the 
improved event selection methods and better-tuned 
MC simulation. 

• Small clustered hits cut 

Using DT data and MC simulation at the position 
(x=-12 m, y=0 m, z—0 m), a difference in cut ef- 
ficiency between DT data and MC events of 2% is 
obtained. Considering the entire fiducial volume 
and the energy region of the solar neutrino signal, 
this difference corresponds to a 0.5% uncertainty 
on the total flux in the 5.0-20 MeV region. A value 
of 2% is assigned to the spectral shape uncertainty 
for the 5.0 to 6.5 MeV bins. 

• Signal extraction method 

The solar neutrino flux is obtained by fitting a solar 
angle distribution (see Figure [25]). To check biases 



Source 


Total Flux 


Energy scale 


±1.4 


Energy resolution 


±0.2 


8 B spectrum 


±0.2 


Trigger efficiency 


±0.5 


Angular resolution 


±0.67 


Fiducial volume (vertex shift) 


±0.54 


Event quality cuts 




- Quality cut 


±0.4 


- Hit pattern cut 


±0.25 


- Second vertex 


±0.45 


Spallation 


±0.2 


External event cut 


±0.25 


Small cluster hits cut 


±0.5 


Background shape 


±0.1 


Signal extraction method 


±0.7 


Livetime 


±0.1 


Cross section 


±0.5 


Total 


±2.1 



TABLE IV: Summary of the systematic uncertainty of the 
total flux in %. 

of the fitting method to the flux value, the solar an- 
gle fit is applied to dummy data which have known 
numbers of signal and background events. As a re- 
sult, ±0.7% difference is found between the input 
and the output number of signal events for the to- 
tal flux, and ±2% difference is found especially in 
the 5.0 < E < 5.5 MeV region. 

1. Summary of systematic uncertainties 

• Flux, time variation, day-night asymmetry 

The systematic uncertainties on total flux, time 
variation and day-night asymmetry are summa- 
rized in Table IIVI The systematic uncertainty on 
the total flux is estimated to be 2.1%. This is 
about two thirds of the corresponding SK-I value. 
The main contributions to the improvement are the 
vertex shift, angular resolution, and event selection 
uncertainties, which are re-estimated for SK-III. 

• Spectrum 

The systematic uncertainty of the spectral shape 
consists of two components: 

— Energy-correlated: 

The energy-correlated systematic uncertain- 
ties are obtained by counting the number of 
events in the solar neutrino MC simulation 
with artificially shifted energy scale, energy 
resolution and 8 B v energy spectrum. The 
results of the calculations are shown in Figure 
[251 These correlations are taken into account 
in the oscillation analysis. 

— Energy-uncorrelated: 

The energy uncorrelated spectral uncertain- 
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FIG. 25: Energy-correlated systematic uncertainties. The 
solid, dotted, and dashed lines show the uncertainties of the 
8 B spectrum, the energy scale, and the energy resolution, re- 
spectively. 



Encrgy(McV) 


5-5.5 


5.5-6 


6-6.5 


6.5-7 


7-7.5 


7.5-20 


Trig cff 


±2.4 


±0.9 


±0.1 








(a) 


±2. 


±1.75 


±1.5 


±1.25 


±1.0 




(b) 








±0.25 


±0.25 


±0.25 


Small cluster hits cut 


±2. 


±2. 


±2. 








External event cut 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


Fiducial volume (vertex shift) 


±0.5 


±0.1 


±0.1 


±0.1 


±0.1 


±0.1 


BG shape 


±0.2 


±0.8 


±0.2 


±0.2 


±0.2 


±0.2 


Sig.Ext. 


±2.1 


±0.7 


±0.7 


±0.7 


±0.7 


±0.7 


Cross section 


±0.2 


±0.2 


±0.2 


±0.2 


±0.2 


±0.2 


Total 


±4.3 


±3.0 


±2.6 


±1.5 


±1.3 


±0.8 



(2.32±0.04(stat.)±0.05(sys.)) x 10 6 cm^sec -1 . 

This result is consistent with SK-I (2.38 ± 0.02 (stat.) ± 
0.08(sys.) x 10 6 cm-^ec- 1 ) and SK-II (2.41 ± 
Q.05(stat.)lg;ig(sys.) x 10 6 cm^sec" 1 ). The SK-I and 
II values are recalculated using the 8 B spectrum of [lOj. 

Figure [2U shows the angular distribution of extracted 
solar neutrino events. 
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FIG. 26: The angular distribution of the solar neutrino final 
sample events. The dotted line seen under the peak in the 
solar direction represents background contributions. 



TABLE V: Energy-uncorrelated systematic uncertainty on 
the observed spectrum shape in %. 



ties are listed in Table [V] To obtain the effect 
on the energy spectral shape, relative differ- 
ences of each uncertainty between each energy 
bin are studied. For example, the uncertainty 
of the fiducial volume due to the vertex shift 
is 0.54% in total, and the relative differences 
are obtained as 0.1% for 5.5-20 MeV region, 
but 0.5% for 5.0-5.5 MeV region with the tight 
fiducial volume. These uncertainties are taken 
into account in the oscillation analysis without 
correlations. 



I. Energy Spectrum 

The recoil electron energy spectrum is obtained from 
5.0 to 20.0 MeV in 21 bins. The definition of the energy 
bins is given in Table VI, which shows the observed and 
expected event rates. Figure [37] shows the observed en- 
ergy spectrum divided by the SSM(BP2004) Q without 
oscillation. The line in the Figure [57] represents the total 
SK-III average (flat data/SSM prediction without neu- 
trino oscillation effect). The \ 2 value is 27.1/20 dof for 
the flat prediction and 26.8/20 dof for the prediction with 
the best fit neutrino oscillation parameters obtained by 
the global solar analysis (see Section HVC[) . This result 
indicates no significant spectral distortion. 



H. Total Flux Result 

Recoil electrons from elastic solar ^-electron scatter- 
ing are strongly forward-biased. SK-III statistically sep- 
arates solar i/'s from background with an unbinned like- 
lihood fit to the directional distribution with respect to 
the Sun. For a live time of 548 days of SK-III data, from 
5.0 to 20.0 MeV, the extracted number of signal events is 
8132±^|f (stat.) ±186(sys.). The corresponding 8 B flux 
is obtained using the 8 B spectrum of 10] to be: 



IV. OSCILLATION ANALYSIS WITH SK-III 
RESULT 

Oscillations of solar neutrinos have been studied by nu- 
merous experiments. The results from such experiments 
have placed increasingly stringent constraints on the mix- 
ing angle between neutrino mass and flavor eigenstates 
as well as on the neutrino mass difference. 

In section IIVB1 the result of the two-flavor neutrino 
oscillation analysis using SK-I, II, III data is presented. 
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TABLE VI: SK-III observed energy spectra expressed in units of event/kton/year. The errors in the observed rates are statistical 
only. The expected rates neglecting oscillation are for the BP2004 SSM flux values. 8 Z is the angle between the z-axis of the 
detector and the vector from the Sun to the detector. 
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FIG. 27: Ratio of observed and expected energy spectra. The 
dashed line represents the SK-III average. 



A conventional two flavor analysis is done in order to 
compare directly with previous results. The results from 
all other solar neutrino experiments and KamLAND are 
combined in section llV CI The two-flavor analysis in SK- 
III is accomplished in basically the same way as the pre- 
vious SK-I and SK-II analyses. Updates to the experi- 
mental data which are used in this analysis are explained 
in each section. 

In section llVDl the results of the three flavor analysis 



are shown. The experimental data have been improved 
statistically and systematically and now have a few per- 
cent uncertainty. In this paper, we present the first result 
of a three-flavor analysis using full SK data information. 

A. x 2 definition 

The oscillation analysis of SK uses the spectrum, time 
variation (zenith angle dependence), and total flux in de- 
termination of the solar neutrino oscillation parameters 
(6*i2, Ami2). For each set of oscillation parameters, the 
total 8 B and hep neutrino fluxes are fit to the data. The 
entire SK-III observed spectrum is utilized for a 5.0 MeV 
threshold. The MSW [15( v e survival probabilities are 
numerically calculated from the solar matter distribu- 
tion provided by SSM(BP2004). The absolute 8 B and 
hep neutrino flux predictions of the SSM are only used 
as a normalization. The uncertainty of the 8 B neutrino 
flux of the SSM is not used in the calculation of the % 2 . 
We use the uncertainty of the SNO neutral current (NC) 
measurement instead (see the next section); hence, the 
absolute 8 B flux predicted by the SSM does not affect 
the fitted 8 B neutrino flux value. For the hep neutrino 
flux, we use the uncertainty of the SSM, since experimen- 
tal uncertainty is still large; hence the fitted hep neutrino 
flux is constrained by the SSM prediction. The predicted 
neutrino spectrum is then converted to an expected SK- 
III rate spectrum by using the v — e elastic scattering 
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cross section and the SK-III detector energy resolution. 
To account for the systematic uncertainties in energy res- 
olution as well as the energy scale and the 8 B neutrino 
spectrum model shape, the combined rate predictions are 
modified by energy shape factors, f{Ei,5B,5s,5n). The 
quantities 5b, <5s, and Sr represent uncertainty in the 
8 B neutrino spectrum, SK-III energy scale, and SK-III 
energy resolution respectively. The function / serves to 
shift the rate predictions corresponding to a given un- 
certainty S in the data rate. The following equation 
describes the SK-III spectrum \ 2 along with energy- 
correlated systematic error shape factors applied to the 
expected rate: 



Xsk-iii 
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(di - (pbi + rihi) x f(Ei,6 B ,5 s ,6 R ))< 




1 10 10' 
tan 2 (0) 
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'SNO 
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(4.1) 



FIG. 28: 95% C.L. allowed region from SK-I,II,III combined 
analysis. The 8 B flux is constrained by the SNO NC rate 
(LETA and phase-Ill). 



where di is the observed rate divided by the expected, 
unoscillated rate for the i th energy bin. Similarly, bi and 
hi are the predicted MSW oscillated rates divided by the 
unoscillated rate for 8 B and hep neutrinos respectively. 
P (rj) scales the 8 B (hep) neutrino flux. L is the unbinned 
time- variation likelihood for the SK-III solar zenith angle 
flux variation above a 5.0 MeV threshold. This likelihood 
is analogous to the one used in SK-I and SK-II. The last 
two terms in the Equation 14 . 1 1 are the constraining terms 
of 8 B and hep fluxes respectively. The value of <^sno 
is the NC value coming from the measurement of SNO 
divided by the SSM flux value, and ctsno is the accompa- 
nying uncertainty of their measurements [19, 20]. The 
numerical values are <^sno = 0.899 and ctsno = 0.032. 
CThcp is the uncertainty of SSM prediction on the hep 
neutrino flux, which is 16% [7|. 



B. Oscillation Results - SK with constrained flux 

The oscillation analysis is performed by including 
X 2 terms corresponding to the SK-I and SK-II values 
(namely, the spectrum and unbinned time variation for 
SK-I and SK-II). By constraining the 8 B flux to the to- 
tal NC flux value from SNO, allowed parameter regions 
can be obtained. Figure [28] shows allowed regions at 95% 
confidence level. The result is consistent with previous 
SK-I and II results. This result is the first to show that 
the energy spectrum and the time variation of the solar 
neutrino flux measured by SK favor only LMA solution 
at 95% C.L., by constraining the 8 B neutrino flux to the 
SNO NC flux and the hep neutrino flux to the SSM pre- 
diction. 



C. Combined Oscillation Results from several 
experiments 



The combination of other solar neutrino experimental 
results such as the SNO, Borexino and radiochemical re- 
sults with the SK combined analysis is accomplished with 
a two-flavor neutrino oscillation analysis by constructing 
a global x 2 ■ F° r the SNO results, the total CC rates ob- 
served in the 306-day pure DoO phase (SNO-I) [13, 391- 
day salt phases (SNO-II) [Hf , and 385-day NCD phase 
(SNO-III) pi, the combined NC rates of LETA [U 
and SNO-III, and the predicted day-night asymmetry for 
SNO-I and II are used. In this analysis, the correlations 
between the SNO CC and NC rates are not taken into 
account. The 7 Be solar neutrino flux of Borexino's 192- 
day 12111 and the radiochemical experiments of Homes- 
take GALLEX-GNO [H, and SAGE [H are then 
included into the global x 2 with the fluxes and their cor- 
relations calculated by SSM in a way shown by (2f| ■ Fig- 
ure [25] shows the combined solar allowed region. The 
best fit parameter set is sin 2 9i 2 = 0.30l ; 01 (tan 2 9i 2 = 
0.42+°;°^) and Am^ = 6.2±^| x 10~ 5 eV 2 , consistent 
with the SK-I global analysis. In addition, combining 
the above and KamLAND data, the best fit parameter 
set is sin 2 6 12 = 0.31 ± 0.01 (tan 2 6 12 = 0.44 ± 0.03) and 
Amlj = 7.6 ± 0.2 x 10~ 5 eV 2 as shown in Figure H3 

Figure [31] shows the 95% allowed region for all solar ex- 
periments before and after the SK-III result is included. 
As shown in the figure, the SK-III result contributes 
about 5% improvement to the uncertainty of Am 2 2 . 
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FIG. 29: Allowed region for all solar experiments at 95% C.L.. 
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FIG. 30: Allowed region for all solar experiments and Kam- 
LAND for two-flavor analysis at 95% C.L.. 



D. Three flavor analysis 

In a three-flavor analysis, the calculation of oscillation 
probability is based on 2l\. The probability can be cal- 
culated with three parameters: 6\2, #13, and AiHj 2 , as- 
suming Am 2 2 <C A77123 ~ Amij. We fixed Am| 3 = 
2.4 x 1CP 3 eV 2 and the normal hierarchy is assumed. For 
the solar neutrino oscillation, the other mixing parame- 
ters are irrelevant, but we set #23 = t/4 and 5c p = in 
our calculation. 

As done for the two-flavor analysis, the oscillation 
probabilities depending on different zenith angles are cal- 
culated, and then the rate of the radiochemical experi- 
ments, Borexino ( 7 Be neutrino flux), and SNO (CC for 



FIG. 31: Allowed region for all solar experiments for two- 
flavor analysis at 95% C.L.. The light green contour and the 
yellow star show the result without SK-III, and the light blue 
contour and the black star show the result with SK-III, which 
is same as Figure [29] 



all phases, NC for LETA and phase-Ill) are calculated. 
SNO day-night asymmetries for phase I and II are also 
predicted. 

The KamLAND spectra and rates from surrounding re- 
actors are calculated based on the information from the 
KamLAND official database [28| and the published pa- 
per [2(| . Because the information in (2|| is used for their 
second result [2{| , we normalized the calculated neutrino 
spectrum without oscillation to that of Figure 1 of (2(| ■ 
The oscillation probability is calculated using three- flavor 
vacuum oscillation. The amount of background and sys- 
tematic uncertainties for each energy bin are read from 
Figure 1 of [26(. The background is fixed in our calcula- 
tion of x 2 - The systematic uncertainty for each energy 
bin is treated as an energy-uncorrelated systematic un- 
certainty. The systematic uncertainty of event rate on 
Table 1 of 20\ is also taken into account. We checked 
first that the contour of the two-flavor analysis in Figure 
2 of [26| was reproduced with 0\% — 0. We also checked 
that our three-flavor contour was very close to the con- 
tour presented in [30|, and consistent with their latest 
result (HI. 

The oscillation parameters are scanned in the fol- 
lowing regions: 10 -5 eV 2 < Am 2 2 < 2 x 10~ 4 eV 2 , 
0.1 < tan 2 6>i2 < 1.0, and < sin 2 13 < 0.25. Figure 
[32] shows the allowed region of the solar neutrino param- 
eters, ((9i2,Am 2 2 ), obtained by the three-flavor analysis 
of the global solar results. The allowed region for Kam- 
LAND obtained by our three-flavor analysis and the al- 
lowed region for global solar and KamLAND combined 
analysis are also shown in Figure[32l Inclusion of another 
oscillation parameter 613 results in a weaker constraint 
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on the solar parameter space. 

Figure I3"3l shows the allowed region in (#12, #13) space 
obtained from the global solar analysis and our Kam- 
LAND analysis. As shown in the figure, in the global 
solar contour, the larger value of #13 prefers the larger 
value of 9i2 , while in the KamLAND contour the larger 
value of #13 prefers the smaller value of 0i2- The 
global solar analysis finds that the best fit values at 
sin 2 0i2 = 0.31±0.03 (tan 2 12 = 0.44±0.06) and Am^ = 
6.0t 2 ; 2 x l(T 5 eV 2 . Combined with the KamLAND re- 
sult, the best-fit oscillation parameters are found to be 
sin 2 6»i2 = 0.31±{J:o2 (tan 2 12 = 0.44™) and Am^ = 
7.7±0.3 x l(T 5 eV 2 . The best fit value of sin 2 9 13 is 0.01, 
and an upper bound is obtained, sin 2 #13 < 0.060 at the 
95% C.L., for the global solar analysis. Combining with 
the KamLAND contour, the best fit value of sin 2 6*13 is 
0.025±g;S| and the 95% C.L. upper limit of the sin 2 9 13 
is found to be 0.059. 




0.1 0.2 0.3 0.4 0.5 

sin 2 «=y 



FIG. 32: Allowed region in solar parameter space (612, Am 2 ) 
obtained by the three-flavor analysis. The thick lines and 
the star mark show the allowed regions and the best fit point 
of the global solar analysis. The thin lines and the square 
mark show the allowed regions and the best fit point of our 
KamLAND analysis. The filled areas and the filled circle 
mark show the allowed regions and the best fit point of the 
combined analysis. For all regions, the innermost area (red), 
the middle area (green) and the outermost area (blue) show 
68.3, 95, 99.7 % C.L. respectively. 

The flux value of 8 B neutrinos can be extracted using 
the oscillation parameters obtained from the fitting of 
the global solar and KamLAND result. As in Equation 
14.11 j3 is a free parameter to minimize the \ 2 an d there 
is no constraint from the SSM prediction in X^sk+sno- 
Table IVIII summarizes the scaled 8 B flux values by us- 
ing j3 m at the best fit point obtained by the global solar 
analysis and the global solar + KamLAND analyses in 
both two and three flavor analyses. The size of the error 
corresponds to the maximum and minimum flux values 
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FIG. 33: Allowed region in (#12,013) space obtained by the 
three-flavor analysis. The definitions of marks and lines are 
same as in Figure [32] 



a B flux (xKfcm-^s- 1 ) 
Global solar (2 flavor) 5.3 ±0.2 
Global solar + KamLAND (2 flavor) 5.1 ±0.1 
Global solar (3 flavor) 5.3 ±0.2 
Global solar + KamLAND (3 flavor) 5.3 ±02 



TABLE VII: 8 B neutrino flux obtained from the oscillation 
parameter fitting. 



among the la oscillation parameter region. As shown in 
the table, the 8 B flux agrees well with the latest SSM 
prediction (32|, and the size of the uncertainty is 2 ~ 3% 
which is consistent with the SNO result [20j |. 



V. CONCLUSION 

Super-Kamiokande has measured the solar 8 B flux to 
be (2.32±0.04(stat.)±0.05(sys.)) x 10 6 cm^sec" 1 dur- 
ing its third phase; the systematic uncertainty is smaller 
than for SK-I. Combining all solar experiments in a 
two flavor fit, the best fit is found to favor the LMA 
region at sin 2 6 12 = 0.30t l (tan 2 6 12 = 0.42±g$) 
and Amfj = 6.2±i;| x 10~ 5 eV 2 . Combined with the 
KamLAND result, the best-fit oscillation parameters are 
found to be sin 2 12 = 0.31 ±0.01 (tan 2 12 = 0.44±0.03) 
and Am|i = 7.6 ±0.2 x 10" 5 eV 2 , in excellent agreement 
with previous solar neutrino oscillation measurements. 
In a three-flavor analysis combining all solar neutrino ex- 
periments and the KamLAND result, the best fit value 
of sin 2 0i3 is found to be 0.025±g;gi| and an upper bound 
is obtained as sin 2 i3 < 0.059 at 95% C. L.. 



19 



VI. ACKNOWLEDGMENTS 

The authors gratefully acknowledge the cooperation of 
the Kamioka Mining and Smelting Company. Super-K 
has been built and operated from funds provided by the 
Japanese Ministry of Education, Culture, Sports, Sci- 
ence and Technology, the U.S. Department of Energy, 



and the U.S. National Science Foundation. This work 
was partially supported by the Research Foundation of 
Korea (BK21 and KNRC), the Korean Ministry of Sci- 
ence and Technology, the National Science Foundation 
of China, and the Spanish Ministry of Science and Inno- 
vation (Grants FPA2009-13697-C04-02 and Consolider- 
Ingenio-2010/CPAN). 



S. Fukuda et al, Nucl. Instr. Meth. A501 (2003) 418. 
J. Hosaka et al, Phys. Rev. D73 112001 (2006). 
J. P. Cravens et al, Phys. Rev. D78 032002 (2008). 
S. Fukuda et al. Nucl.Instr.and Meth. A 51 418 (2003). 
M. Smy, Proc. of the 30th ICRC, 5 1279 (2007). 
R.M.Pope and E.S.Fry, Appl. Opt. 36, 8710-8723 (1997). 
J.N.Bahcall and M.H.Pinsonneault, Phys. Rev. Lett. 92 
121301 (2004). 

C.E.Ortiz,A.Garcia,R.A.Waltz,M.Bhattacharya,A.K.Komives[25 
Phys. Rev. Lett. 85 2909 (2000). 
J.N.Bahcall et al, Phys. Rev. C54 411 (1996). 
W.T.Winter,S.J.Freedman,K.E.Rehm,J.P.Schifrer, Phys. 
Rev. C, 73 025503 (2006). 

J.N.Bahcall et al, Nucl. Phys. B (Proc. Suppl), 77, 64-72 
(1999). 

M.Nakahata, et al. Nucl. Instr. Meth. A421 113 (1999). 
E.Blaufuss, et al. Nucl. Instr. Meth. A458 636 (2001). 
R. M. Pops et al. Applied Optics 35 33 (1997). 
S.P. Mikheyev and A.Y. Smirnov, Sov. J. Nucl. Phys. 42 
913 (1985) ; L. Wolfenstein, Phys Rev. D17 2369 (1978). 
M.B.Smy, arXiv:hep-ex/0202020. 
B. Aharmim et al, Phys. Rev. C, 75 045502 (2007). 



[18 
[19 
[20 
[21 
[22 
[23 

[24 



[26 
[27 

[28 

[29 
[30 

[31 
[32 



S.N.Ahmed et al, 
B. Aharmim et al. 

B. Aharmim et al. 

C. Arpesella et al, 
B.T.Cleveland et 
M.Altmann et al, 



Phys. Rev. C. 72 055502 (2005). 
Phys. Rev. Lett. 101 111301 (2008). 
Phys. Rev. C 81 055504 (2010). 
Phys. Rev. Lett. 101 091302 (2008). 
il, Astrophys. J. 496 505 (1998). 
Phys. Lett. B490 16 (2000); also pre- 
sented at Neutrino 2004 by C.Cattadori. 
J.N.Abdurashitov et al, Phys. Rev. C60 055801 (1999). 
G.L.Fogli,E.Lisi,A.Marrone,D.Montanino,A.Palazzo, 
Phys. Rev. D66 053010 (2002). 
S.Abe et al, Phys. Rev. Lett. 100 221803 (2008). 
V.Barger,K.Whisnant,S.Pakvasa,R. J.N. Phillips, Phys. 
Rev. D. 22 2718 (1980). 
KamLAND official data base wep page: 

http: //www. awa.tohoku.ac.jp/KamLAND/datarclcasc/2ndrcsult. html 

T.Araki et al, Phys. Rev. Lett. 94 081801 (2005) 
K.Ichimura for KamLAND collaboration, 34th Interna- 
tional Conference on High Energy Physics (ICHEP 2008). 
A.Gando et al. II arXiv: 1009.4771 1 
A.M.Serenelli et al. arXiv:0909.2668 v2 



